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© The screen is provided by surrounding the coil producing the magnetic field with a set of electrical 
conductors the currents within the conductors being controlled in such a manner that the field is neutralised in a 
specific region of space, the current distribution within the conductors being determined by calculating the 
current within a hypothetical superconductive shield which would have the effect of neutralising the field, the 
current through the conductors thereby- being a substitute for the superconductive shield. 
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MAGNETIC FIELD SCREENS 



This invention relates to magnetic field screens and has application in NMR imaging apparatus. 

Current carrying magnet coils are used for a variety of purposes in NMR imaging apparatus. Examples 
include large electro-magnets designed to provide static magnetic fields to polarise nuclear spins, magnetic 
field gradient coils which superimpose gradients onto the static polarising field and rf transmitter and 
5 receiver coils. 

. In many cases the design of a magnet coil is such as to optimise the magnetic field within a desired 
volume. However the coil inevitably produces an extraneous magnetic field outside that volume, especially 
relatively close to the coil. In the case of large bore static electromagnets the high fields they generate will 
produce undesirably strong extraneous fields at distances outside the magnet corresponding to many 

w magnet diameters. Such magnet systems therefore require much free and unusable space around their 
installation. Stray magnetic fields may be intolerable in hospitals because of iron structures positioned 
around the installation site which vitiate the homogeneity of the magnetic field Additionally, electronic 
equipment may not perform well in an environment which has an extraneous magnetic field. 

Furthermore, most NMR imaging systems utilise rapidly switched magnetic field gradients in their 

/s operation. A major problem epecially where super-conductive magnets are used, is the interaction of the 
gradient field with the magnet itself. Existing attempts to minimise this interaction include the use of 
conducting metal screening sleeves. However, induced currents in these sleeves or in the heat shield of the 
magnet decay with uncontrolled relaxation times which make it difficult or even impossible to implement 
some of the faster and more efficient NMR imaging methods. This is because the decaying current 

20 produces image fields superimposed on the desired gradient field. This uncontrolled time dependence 
introduces phase artefacts which can completely ruin the image. 

In order to provide adequate access for patients, and to improve gradient uniformity, it is desirable to 
maximise the diameter of the magnetic field gradient coils in an NMR imaging machine. However, this often 
causes the coils to be close to other conductors, either the surfaces of cryogenic vessels (in superconduc- 

25 ting magnet systems), electromagnetic coil supports (in resistive magnet systems), or ferromagnetic pole 
pieces (in ferromagnetic systems). When gradients are switched rapidly, as with many imaging techniques, 
eddy currents are induced in these conductors which then contribute additional field gradients varying in 
time and potentially vary non-uniform in space. Typical time constants for the decay of the eddy currents 
vary from a few milliseconds to hundreds of milliseconds, depending on the type of main magnet and the 

30 specific coil configuration. 

The commonest solution to this problem is to tailor the input applied to the amplifiers generating the 
gradient coil currents in such a way that the gradient fields themselves follow the prescribed time variation. 
The input voltage and gradient coil currents are characteristically over-driven for the initial part of the on- 
period of the gradient. But this remedy has a major disadvantage. If the gradient coils are placed close to 

35 the coupled surfaces, so that the eddy current field gradients may have the same uniformity as the desired 
gradient the gradient coils become very inefficient and a large over capacity in the gradient current 
amplifiers is required, since the 'reflected* fields will be targe and in the opposite sense from the desired 
fields. If, on the other hand, the gradient coils are reduced in size, to avoid the amplifier capacity problem, 
then the reflected gradient fields will in general be non-linear over the region of interest. Furthermore, in 

40 either case there are likely to be reflected fields from more distant conductors in the main magnet structure, 
each with its distinct time constant and spatial variation. 

The only effective solution is in some way to reduce the gradient fields to zero at a finite radius outside 
the coils, so that no undesired conducting structures may be coupled to them. 

Partially effective methods for magnetic screening in specific coil geometries have been proposed 

4$ hitherto in particular U.S. Patent Nos. 3466499 & 3671902. These geometries are not generally useful in 
NMR and NMR imaging. 

It is an object of the invention to provide more efficient and effective magnetic field screens for coil 
geometries useful in NMR and NMR imaging. 

It is a further object of the invention to provide efficient magnetic field screens for any coil design, 
so It is a still further object of the present invention to provide magnetic field screens for any desired 
component or components of a magnetic field. 

According to the invention a screen for a magnetic field comprising a set of electrical conductors and 
means for supplying the conductors of the set with electrical currents of magnitude such that a) the 
resultant current distribution approximates to the induced current distribution in a hypothetical continuous 
superconductive metal surface positioned in the place of said set so as to appear as a complete reflector of 
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magnetic field, and b) the resultant current distribution in this or other screens behaves alone or in a 
combination with said other screens in such a way as to appear to selectively reflect and/or transmit desired 
components of magnetic fields of specific configuration through said screen or screens. 

Preferably the current distribution localised to the surface of a hypothetical conducting sheet or sheets 
5 is determined by the deconvolution of the magnetic field response function of the unit line elements of that 
current with the field to be screened; such problems being most conveniently solved in reciprocal space, 
which is defined by those co-ordinates conjugate to real space used in appropriate integral transform. 

More preferably the problems are solved in Fourier space which is a particular example of reciprocal 
space. 

io The present invention also provides a method of designing a screening coil for selectively screening the 
field of a magnetic coil. 

The present invention further provides a gradient coil system for use in NMR apparatus including a 
main coil designed to provide a gradient field and a screen coil surrounding the main coll. 

In one preferred arrangement the conductors of the set are regularly spaced apart from each other. 
is They may be connected electrically in parallel and have different values of resistance in order to produce 
the desired current distribution. In embodiments of the invention the different values of resistance of the 
conductors may be produced by different thicknesses of the respective conductors or constructing them 
with different compositions having appropriate values of electrical resistivity. 

In alternative preferred arrangements the conductors of the set carry equal currents but are spaced 
20 apart from each other by different spadngs so as to produce the desired current distribution. 

It is a further object of the present invention to reduce acoustic vibration in MR gradient coils by using 
the active screening hereinbefore described. 

In examples of carrying out the above invention a set of conductors are arranged on a cylindrical former 
and with appropriate spacing can be fed with equal currents. The spacing of the conductors of the set and 
25 the magnitude of the currents are calculated using the currents induced in a flat superconductive screen as 
the starting point. 

Additionally according to this invention the theoretical current distribution in a continuous superconduc- 
tive cylindrical shield positioned in the place of the aforementioned set of electrical conductors is calculated. 
Such a calculation enables an improved screen to be provided, especially where a cylindrical screen is 
30 required. 

The calculations described herein represents an analytical formulation of the problem, enabling a fully 
general calculation of the current density in a cylinder required to cancel outside the cylinder the magnetic 
fields generated by coils inside. The results obtained are firstly applied to passive shielding, using a thick 
high-conductivity or superconductive cylindrical tube to solve the reflected fields problem without sacrificing 
35 gradient uniformity. The conclusions arrived at can then be applied to active shielding; the calculated 
current densities within the skin depth of a thick cylinder or in a superconducting cylinder are mimicked 
using a suitable arrangement of a set of wires supplied with currents of appropriate magnitude. 

For passive shielding using a normal conducting cylinder to be effective, the skin depth a in the shield 
at the lowest frequency represented in the particular gradient switching sequence must be much smaller 
40 than the thickness d of the wall of the cylinder. For an echo planar switched gradient, for instance, with a 
fundamental frequency of 1 KHz, this entails a wall thickness of-10 mm of aluminium. For switched 
gradients such as that used for slice selection, where there is a non-zero d.c. component of the fields at the 
cylinder surface, passive shielding is not appropriate. 

When the criterion of b /d <*: 1 is met. the time dependence of reflected fields will be identical to that of 
4$ the applied field, and it is only the spatial non-uniformity of the net magnetic field which is of concern. 
Analytic calculation or the ensuing fields enables the necessary corrections to coil spacings to be made. 

By Lenz's law, a magnetic field screen constructed in accordance with the above design criteria 
produces a magnetic field which opposes the field generated by the primary magnetic coil that it is 
designed to screen. For a given current in the primary coil the resultant magnetic field generated within the 
so volume embraced by the coil is reduced and its spatial variation is also changed by the presence of the 
screen currents thus introducing undesirable variations in the primary field. 

It is therefore a further object of the invention to provide a screening coil arrangement in which the 
above disadvantages are overcome. 

Accordingly the present invention also provides a magnetic coil surrounded by two active magnetic 
55 screening coils, namely an inner screen and an outer screen, each respective screen comprising a set of 
electrical conductors and means for supplying the conductors of the set with electricaJ currents of 
magnitudes such that there is no appreciable magnetic field outside the outer screen and the field within the 
inner screen substantially corresponds to the field that would be provided by the said magnetic coil if the 
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screens were not present 

In order that the invention may be more fully understood reference will now be made to the 
accompanying drawings in which:- 

Figures 1 to 10 are explanatory of the underlying theory; 
5 Figures 11 to 17 illustrate in diagrammatic form various embodiments of the invention; 

Figure 18 shows a cylindrical co-ordinate system used in calculating the magnetic fields produced by 
current flow on a cylindrical surface; 

Figure 19 illustrates saddle-shaped coils used as magnetic gradient field coils; 
Figure 20 illustrates the configuration of one octant of a set of screening coils in accordance with the 
io calculations described herein: 

Figure 21 shows curves for optimising the positions of the arcs of the saddle coils of Figure 19; 
Figure 22 illustrates diagrammatically the magnetic field produced by an unscreened primary coil in 
the form of a singfe hoop; 

Figure 23 illustrates the same hoop with a double screen embodying the invention; 
75 Figure 24 is a graph showing the magnetic field at different radiai positions that is produced by the 

Figure 23 arrangement; 

Figure 25 is a perspective view of a double saddle coil used to produce transverse gradient fields 
which is screened by a double screen embodying the invention: 

Figures 26 to 33 illustrate ways of providing varying current distribution over a desired area; 
20 Figure 34 shows a parallel band arrangement; 

Figure 35 shows a series arrangement; 

Figure 36 shows various arrangements for assisting in reduction of acoustic vibration in MR coils by a 
gradient wire with (a) a wire pair (b) a wire array (c) a pair of conducting plates used to nullify the main field 
in the neighbourhood of the gradient coil and (d) a double active screen arrangement; 
25 Figure 37 shows winding strategies for screening wires about gradient coils of the (a) standard 

circular design and (b) saddle design; and 

Figure 38 illustrates the screening of one coll from the magnetic field of another. 
The basic theory and ideas are developed from the simple case of an infinite straight wire parallel to an 
infinite flat conducting sheet. Figure 1 shows a long straight wire carrying a current I. The magnetic field B 0 
. ao at a point P which is at a distance r normai to the wire is given by 



35 



a *V 21 

o l 4i' r (1) 



If the current is changing with angular frequency « and the wire is near to an infinite conducting sheet as 
shown in Rgure 2 (also if the current is static and the sheet is superconducting), the magnetic field 
undergoes distortion at the metal surface. Let us assume no field penetration into the screen, i.e. the 
40 angular frequency « and the metal conductivity a are sufficiently high, then the boundary conditions for the 
magnetic field at the surface are 
B* * 0 

B, - 0 (2) 

^ The details of the field at a point P due to a wire at a distance d from the sheet may be conveniently 
calculated as shown in Figure 3 using the method of images which assumes a wire at a distance d from the 
other side of the sheet carrying a current -I. 

In general, current in the sheet surface is directly related to the tangential field H y . The total field B 0 at 
P, Rgure 4. is given by Equation (1) and may be resolved into components B* and By in the x and y 
directions respectively 
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3 - & 2L«so (3) 



Taking the image current -I into account, Figure 5, we obtain the total y direction field Btot at point P for a 
/s current I at point 0 and its image current -I at a point 0\ 



- image 



(5a) 
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(47) ^ 1 cos 




(5b) 



We also note from Figure 5 that 
y = r sine* = r' sinfl (6) 
and r cos* + r cosfl = 2d (7) 

If P is on the conducting sheet surface r « r' in which case, see Figure 6. 

B 7 ~ 18) 



1* ■ d 2 + y 2 (9) 
and cos* = d/r (10) 
• From Equations (9) and (10) we obtain 



Consider now the line integral of the magnetic field in and close to the metal surface, Figure 7. By 
Ampere's theorem we have 
0H y .dl = /J y .ds (12) 

where" J y is~the"surface current density and ds an area element and dl a path element For a short path I 
^ parallel to the surface H y constant The line integral is therefore 
Hyl-H v l + (H x + H M )dl * J y l (13) 
But in the metal H/ ■ 0 and H x = H x yielding for dl-K) 

Hy = Jy (H) 
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The surface current distribution is therefore 

aid 



<TT> 



J t (<i) 



(15) 



This function is plotted in Figure 8. 

To determine the field produced by the surface current density distribution, let us assume we have a 
surface distribution Jy within a flat metaJ sheet as in Equation (15). Consider an element of surface dl with 
current 5i. Figure 9. This current is given by 
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rdy 



(16) 



'5 The elemental field at point P which is a distance r from the element and distance d from the sheet is 

li (17) 



6H 



1st > 



20 The tangential component of which is> 

£i sLa6 <5i d 



25 



yd « 



ixr 



(18) 
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2l&~, 1. dv 



(TT> 



1 4</( d 2 * y 2 ) Z ' 



(19) 



The total field is 



40 



which yields 



a y * 2 



2 Id 2 J 



dv 



w o 



(tt) — 2 t 2 



(20) 



(21) 



Thus the field at point P distance d from the surface current distribution is equivalent to a mirror current of 
•I at distance 2d from P. 
50 The totaJ surface current 



I y - I J..dy = 1 <22) 
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The results derived above suggest that instead of using a metal plate to screen alternating fields, an 
active screen comprising a mesh of wires may be used in which a current pattern is externally generating to 
mimic a desired surface distribution. This situation is shown in Figure 10. Wires at positions yi, y 2 ....y n in 
the y direction alt carrying an equal current i are spaced at different discrete intervals corresponding to 
5 equal areas A under the J t curve of Figure 8. For equal currents in the wires, the wires must be unequally 
spaced with positons y„ which satisfy the relationship 



(«•»>- ?fw*y (23a , 

for an even array of 2N wires spread about the y origin. For an odd array of 2N + I wires with one wire at 
J5 the origin we have 

20 0 



25 



30 



where n ■ 1, 2, 3... 

Alternatively, wires may be equally spaced with different currents chosen such that 

^(a^y)dy (24) 

(n-1) 



In either case, since we have arranged to satisfy the original boundary conditions on the surface of a 
fictitious plate, the magnetic field in the half plane (x.ty) approaches zero. The degree of screening 
depends ultimately on the number of wires used in the active screening mesh. An example of such an 
active screen is shown in Figure 11 for a distribution of current carrying conductors corresponding to Figure 

35 10 all carrying an equal current L The field B on the opposite side of the screen is effectively zero. 

If two parallel screens are used, each will have primary current distributions of Ji(di) and Ji(d2> given 
by Equation (15). This is shown in Figure 12. However, each induced current distribution will induce further 
changes in the distribution in the opposite plate represented by additional terms J n (d n ). This is equivalent to 
a set of multiple reflections. Figure 13, which correspond to an infinite set of images. When di = 62 s d, 

40 image currents occur at x ■ i2nd, n ■ 1, 2... The total induced surface current in each sheet is the sum of 
J(d) *£„J„(d n ). 

The boundary conditions at the metal surface ensure that the normal laws of reflection apply. However, 
when active current screens are used, the reflection laws may be selectively changed to reflect any images 
using the J„(d)'s corresponding to particular distances d. 

45 An example of a gradient active screen is shown in Figure 14 for a single circular hoop of diameter 2a. 
Let this be screened by an active current mesh in the form of a cylinder of radius a + d. in a metal 
cylinder, image hoops appear at radii r = b + d, 3b±d, 5b±d, etc. However, since the effect of these distant 
images diminishes quite rapidly, it is reasonable to approximate a screen with Ji(d) corresponding to the 
plane sheet case, Equation (15). Better approximations may be obtained by an iterative numerical approach. 

50 Although exact solutions for the surface current in a cylinder exist, when actual wire screens are 
constructed, the numerical approach is preferred since it automatically takes account of the finite number of 
wires and their discrete spacing. The screen current is 
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where N 2 is the number of screen wires in the mesh of radius r 2 each carrying a current and Ni the 
number of turns in the primary coil radius n carrying a current li . The factora is of the order of unity and is 
chosen to optimise the screening. The whole optimisation procedure is accomplished by a computational 
process which generalises to applying minimisation of the total field over a limited region of space. 

5 Mathematically this is conveniently achieved by a least squares method. For practical ease, it is desirable to 
have both coil and screen driven from the same current source. Since the total screen current is less than 
the primary current it will in general be necessary to take parallel combinations of the screen mesh so that 
the total system may be driven in series. However, a parallel arrangement is also possible in which the 
screen wires are varied in resistance and/or impedance in such a way that when driven from a voltage 

to source appropriate currents flow. Because of the screening effect the inductance of both versions of the coil 
system should be low. 

Sets of screened hoops may be used to construct screened magnets producing uniform magnetic 
fields. The presence of the screen around one hoop is approximately equivalent to an image hoop with 
negative current producing an opposing field. For a Helmholtz coil pair the intercoil spacing equals the coil 

i6 radius a. When screened, however, the spacing must be changed so as to make the second derivative of 
the field B, with respect to z vanish for the combined coil system. Similarly for a Maxwell pair designed to 
produce a linear magnetic field gradient, the intercoil spacing is ideally a 3. This is shown in Figure 15. The 
two hoops forming a Maxwell pair are screened by a concentric pair of screening meshes offset axially from 
each other. The combined screen current distribution is also shown. Again, however, when screened coils 

20 are used, a new spacing obtains which makes the third derivative of B, with respect to z vanish for the total 
coil system. The process of optimisation of coil geometry can be done analytically for simple coil structures 
as discussed above. For more complicated systems such as screened saddle geometry gradient coils, it is 
preferable simple to find by computational means the position where the chosen derivative vanishes, or is 
minimised. 

25 For NMR imaging sytems using superconducting magnet coils it is convenient to use saddle coils to 
produce the transverse gradients 

dB^ anddB . 

z_ z 

30 a x a y 

In some imaging techniques, at least one of the gradients can be very large making interaction with the 
main magnet potentially serious. Figure 16 shows one half of a G* screened gradient coil system. A 

35 screened single saddle coil is shown in end view at (a) and plan view at (b). The dotted lines correspond to 
the screening mesh. To a first approximation, the screen current profile is Ji(d). Better screening may be 
obtained by an iterative procedure which minimises the field outside the screen. 

Figure 17 shows a screened rectangular Gx gradient coil set with 1/2 screen. Again, if a d«:2a, Jt(d) 
may be used as a good approximation for the screen current distribution. For better results other reflections 

^ can be included or the iterative procedure used to minimise fields outside the coil. 

NMR imaging systems require rf coil systems to deliver rf pulses to the specimen and to receive 
signals induced in the sample. Because of the number of other coil systems required for field gradients, 
there is always a problem of space. With normal rf coil arrangements the field outside the coil drops off 
rather slowly. In order to minimise coil interactions which can lower the Q value, rf coil diameters are often 

^ chosen to be around 0.5 to 0.7 of the gradient coil diameter. With shielded rf coil designs, it may be 
possible to utilise more of the available space without loss of performance. 

A systematic procedure for reducing extraneous magnetic fields outside the active volume of static 
magnets, field gradient coil systems and rf coils has been described. In NMR imaging, reduction of stray 
fields in all three types of coil structure is extremely important. The method utilises active magnetic screens 

M and has the advantage that such screens operate independently of frequency down to dc. Some price is 
paid in terms of reduction of field in the active volume compared with that of the free space value. With 
time dependent gradients, the price is in general acceptable since for fast NMR imaging schemes 
employing rapid gradient switching, active coil screening may be the only way in which such imaging 
schemes may be made to operate in the relatively close confines of an electromagnet 

^ In what has previously been described, iterative and least suqares approximation methods are used to 
obtain actual screening wire positions. It is possible to obtain these positions directly using analytical 
methods. If the gradient coils and screens are located on cylindrical formers, it is natural to use cylindrical 
co-ordinates p, z in order to retain the symmetry of the system. The z axis is taken to lie along the axis 
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of the cylinder as shown in Figure 1 8. 

The vector potential A. is used to describe the magnetic field. This has the components Ap, A$, A z 
given by 

t /(-) dv 1 sin(d - o' ; 
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15 . U, . 



J^.Cr 1 } dv' coa(d - ) 



(26) 



(27) 



where J is the current density and dv is a volume element corresponding to the position vector r\ There is 
no current flow in the radial direction in many problems of interest so J has only z and * components. 

It has been assumed that currents induced in the shield are confined to the surface of a cyiinder of 
radius b. The gradient coils to be shielded are mounted on a cylindrical former of radius a which is coaxial 
30 with the shield. The currents can then be written as 
J = F(z, *) 5 ( P - a) + f(z, *) 6 (p - b) (28) 

where F describes the current in the gradient coils, f describes the current induced in the shield and 5is the 
Dirac delta function. It is possible to derive relationships between F and f from the condition that the radial 
component of the magnetic field is zero on the surface of the shield. The other constraint which is used is 
36 the equation of continuity, which, in the absence of charge accumulation is 
v . J = 0 (29) 

The vital step in the analytic treatment of this system is the use of the Green's function expansion 
« jf^.l^llp.M*-*') =•) < (30) 

where m is an integer p< ( p *) is the lesser (greater) of p and p\ and i m (z) and K m (z) are modified Bessel 
^ functions. To use equation (30) in equations (25), (26) and (27) it is helpful to define a type of Fourier 
transform of f and F as follows: 
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: * U ) . ^ ra t e ' iM 7 cs= e _lkZ f.(«. -) (31) 



TO 



20 



25 



30 



35 



r « (!<) ■ e' lmo 7c= e" U = r (e. =) (32) 

where the quantities F^k) and F^k) are defined in an analogous way. The components of A thus 
become (e.g. for p> b): 
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[M (kS)f^(k) - al (ka)r\'k)j 



(33) 
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r «H (ks 'V,(^))^) 



S " 3* J. to" 1 " . lkt [«i (lrl (») V| ( fc p, 



Similar expressions can be obtained when p < a or a <p<b. 
55 Now let the boundary condition that the radiaJ component of the negative field Bp ■ 0 at p = b be 
applied. This is equivalent to 
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0 oO 



0-0 



3-D 



(36) 



10 



15 



Equations (33) and (34) are now used and terms varying as e 1 "*^ are equated. This gives 

iS Km(kb) [bl m (kb)f ? M + al m (ka)F 7 (k)] 

* fTMbtl^Ckb^Ocb) + l m *,(kb)]K m . 1 <kb)] 

+ FSMaCUOcaJK^k^ + l m *,(ka) K„, M (ka)] (37) 

This equation can be simplified using the relations derived from the equations of continuity. 

± il* . _ Hi 

o 3* - is ( 38) 



for the currents in the shield, and 

I iii ill 

20 3 3 * * z (39) 



for the currents in the gradient coils. These equations are equivalent to 

»:« = -#f7w (40) 

25 and 

FTM = - A*FT(k) (41) 

The recurrence relations for Bessel functions are also used, from which can be derived the identity 

Ul(2.) K^te) + U,(Z,) K^fo) + (2m*/z,z 2 ) l m (zt)K m (z 2 ) = -21 m '(z,) Cte) (42) 
Here the prime denotes the derivative. 
30 Equations (37), (40), (41) and (42) can be combined to obtain the elegant expressions, 

a* I '(ka) ____ 

35 * 



_ a I ' ( ka J 

3 

(44) 



^ These identities provide the means of calculating fields due to any combination of currents constrained 
to flow on the surface of a cylinder inside a conducting shield. This will now be illustrated with an example. 

For transverse field gradients it is common to use a saddle coil configuration such as shown in Figure 
19. Two pairs of saddle coils are shown. One pair extends between -d t and -d 2 along the z axis. The other 
part extends between +di and + d 2 along the z axis. With energisation of the coils as marked a gradient 

so fi sW is produced in the form of a magnetic field along the z axis which has a gradient in the x direction. The 
field produced by such coils extends widely outside the cylinder on which the coils are wound unless they 
are shielded by a conductive sleeve. 

For the standard coil geometry with 120° arcs for the saddle coils. 
f 4,(4 *) * I(t(z.di) + 5(z + di) - 5(z-d 2 ) - *(z + d 2 )} {H(* +5 )(1 - H(+ - § )) - H<* - §2 )(1 - H(+ + 1= )) 

55 ^ 

where H(x) is the Heaviside step function. This has the Fourier transform 
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r.(k) - — — (cosCkdi) - cos(kd,}) 

0 mZ 



x (l - e ia ~) (46) 

This is zero for m even or an integer multiple of 3. This leaves non-zero terms for m = 1 . 5, 7, 1 1 etc. 
'0 The position of the screening arcs can now be determined as follows. First it is necessary to specify the 
continuous current distribution to which the arcs approximate as set out in equations (43) and (44). 
The actual surface current density may then be written 



20 



25 



30 



9 a f JL " 3 



I '(ka) 



J coskz dk (coskdj - coskd,) ' ■ 
6 J - (k0) (47) 



and 



1 m.i Fs.7 sln "t^n mo 



? stnk= < k (C03kd. - eogKC.) V (ka) 
o ok I -(ko) 



(48) 
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The next step is to determine the stagnation point of this current density, that Is, the point S = (0,d) at 
which and f^ and f 2 are both zero, around which the surface current flows. By symmetry this occurs at $ s 
0, and since f z (O, z) = 0 for ail z, it may be found by solving f$ (0, z) = 0, by successive approximation. 

Having found S. the integrated surface current l t crossing the line SP between S and an arbitrary point 
P on the cylinder (Rgure 20) is given by 



l t (0,z) - £r (0,-)cz - b J f 2 (*.s)d« 
(using the surface version of the divergence theorem). Hence 
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10 



V*' C) o a r 3ln ft ?1 rcoskd," - coskd,) V ( * a) . 
; a-»h.7 n 6 dK * V Ub) 

(sinkz cosm^S - sinkd) (50) 

The contours of the induced surface current may now be found by setting l t ($, z) = constant These 
may be translated into positions of screening arcs in the following way. The total current in the cylinder 1,(0, 
0) is divided by 2N t where N is the number of screening arcs required, and the M'th arc consists of the 
contour where 

is V°' 0? 

I^o.z) — Uk - ») (51) 

In practice !,($, z) is calculated over a grid of 50 x 45 points, and the contours are found by linear 
20 interpolation between the calculated points. 

To check that these arcs indeed provide adequate screening of saddle coil fringe fields, the total field 
can be calculated using the Biot-Savart Law, taking as line elements dl the intervals between successive 
calculated co-ordinate pairs along each arc. 

Given a set of saddle coils with radius 0.31 m and arc spacings from the centre di and d? of 0.108 m 
25 and 0.404 m, the shielding produced at 0.55 m radius, using six screening arcs on a cylinder of radius b = 
0.45 m, is as follows: 

Maximum unscreened field = 0.6 x 10"* T/A-turn 
Maximum screened field (z < 0.5 m) = 0.2 x 10" 7 T/A-tum 
Maximum screened field (z < 1.0 m) = 0.47 x 10~ 7 T/A-tum 
30 For comparison, the maximum field at 0.55 m radius produced by an unscreened small scale saddle coil 
set (a = 0.16 m. di = 0.56 m. d 2 = 0.206 m) is 0.66 x 10~ 7 T/A-turn. 

Figure 20 shows the configuration of one octant of the set of screening coils calculated above. 

A larger number of screening arcs, or the use of foil rather than wire conductors, will further reduce the 
fringe fields. 

35 Screening a set of saddle coils, the spacing of which has been optimised without the screen present, 
inevitably reduces the uniformity of the gradients produced. The uniformity may be recovered however by 
adjusting the arc spacing as follows: 

The z component of the magnetic field can be derived from equations (34), (35), and (44) to give 
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45 



K 1 (kb) I 1 (ka) 

° (52) 



50 This becomes, with substrtion of equations (45) and (46) 
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is 



20 



5 - " * -,-1 ; T i- <3k k COS 01(6 C03(kz) I (kp) 



«• 5 1 f! (hit/"' ) / \ i 

x 7 » ' 1 l C8a < ka i> " coa(kd,)J|K «(k«) 



K '(ko): '(ka) 



(53) 

It is now possible to optimise the gradient linearity by adjusting the arc positions of the saddle coils. The 
terms for m = 5 are of fifth order in x. y or z. whereas there are terms form « 1 which are of first order in 
x and of third order in x, y and z. 

The optimum choice of the parameters di and d2 is the one which removes the third order terms. This 
gives the condition 

25 I dk k* (costtcd,) - co3(kd a ;)(KV(ka) - K '^u D y (Wa) 1 ' C 



30 (54) 

There are now two unknown quantities, Dt(=di/a) and D2( s d2/a) and only one constraint, so it is not 
possible to give a unique choice of parameters. However it can be ensured that each parameter separately 
3 8 satisfies the equation 

40 

(55) 

(with o= b/a) for this automatically satisfies equation (54). Values of Ot and D 2 as a function of a are shown 
^ in Figure 21. 

It is possible to improve on this arrangement by altering Di and D2 slightly (subject to the constraint 
given by equation (55)) and minimising the fifth order terms. The final choice of values of Di and D 2 
depend on whether the x or 2 variation is considered to be more Important. The values of D\{a ) and D 2 ( a ) 
shown represent an excellent starting point in the search for the optimum position of the saddle coils. 
^ Referring now to Figure 22 there is shown therein a magnetic coil in the form of a single wire hoop I of 
radius a carrying a current +1. In Figure 23 the same wire hoop I is surrounded by two active magnetic 
screens S1 and S2. Each screen comprises a set of electrical cunent carrying conductors but for simplicity 
the screens are shown as sections of cylinders. Outer screen S1 is a cylinder of radius b and inner screen 
S2 is a cylinder of radius c. 

55 With appropriate screen current density distributions screens S1 and S2 act together as a flux guide 
confining the field lines in the manner indicated. The design criteria for the current density distributions are 
set out below. 

For a single conducting screen S, the boundary conditions of the magnetic field B(r-r') at the surface of 
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the screen require only that the axial component B. (r-r', z'z) is considered which for a coaxial hoop is 
independent of azimulhal angle for perfect screening it is required that 
IB^r-a.z) + B* (r-b.z) = 0 (56) 

for r > b and for all z. where B t (r-a,z) is the primary hoop field per unit current and B| (r-b.z) is the totaJ 
s field generated by the screen. The screen field is the convolution of the surface current density j*(z) with 
the hoop magnetic field response per unit current Equation (56) may therefore be written as 



Note that F?(q) = J:(q);F?(q) = J?(q) 
is f ?ta> = iJ(Q):f?(q) » j T (q) where q may be taken to be k or z. 

The spatial Fourier transform of equation (57) gives 

I B^r-a*) = -B z (r-b,k)j<>(k) (58) 

where k is the reciprocal space wave number. 

Equation (58) may be generalised for two screens with current densities jj (k) and jj (k) as set out 
20 below. For zero field in the range r £ b and the unperturbed hoop field for r S c there is obtained 

I B^k) = -B x (r-b,k)ji (k) - B^r-c.k) jj (k) for r * b (59) 

0 = B,(r-b.k)ji (k) + B^r-ck)^ (k) for r S c (60) 

The current hoop fields are evaluated numerically along the z-axis on the appropriate cylindrical surface 
using the Biot-Savart Law, then Fourier transformed to k-space. This allows numerical solution of the 

2S simultaneous equations (59) and (60) yielding the k-space current densities. These are then inversely 
transformed to yield the actual screen current density distributions. 

In this example consider a primary hoop of radius a = 0.5m shielded by two active screens St, S2 with 
radii b * 0.75 m and c = 1.0 m respectively. The hoop current is 1A. Using the computed distributions, the 
total magnetic field generated by the screen system in the hoop plane (z = 0) is calculated as a function of 

■a, r. This is shown in Figure 24 which is a graph of a calculated z-component of magnetic field B z (r,0) against 
radius in the plane of a doubly screened flat current hoop of dimensions as in Figure 23. The unscreened 
hoop field is shown by the broken curve and is equal to B, for r < 0.75 m. As expected for r < c, the field 
B 2 exactly equals that of the unscreened hoop. For r > b. the field B is zero. Between the screens, the field 
is wholly negative. The total screen currents are given by I, = (0) and I 2 = jj (0). In this case l 2 = -h = 

3S 0.57 I. The results have been produced by numerical methods using a computer. 

The inner screen S2 behaves like a complete absorber of the primary field. However, once trapped 
between the screens, the field is completely internally reflected by S1 and S2. The screen S2 behaves like 
a perfect one-way mirror. Practical screens having these properties clearly cannot be continuous metal 
surfaces. Externally driven discrete wire arrays which approximate the calculated continuous surface current 

40 densities are used instead. For a common screen cunent, wires are placed at positions corresponding to 
equal areas under the current density distribution integrated over discrete intervals, as described in the 
above patent specifications. 

Although active screening of a single hoop has been described in the above example, the principles of 
double screening apply to other geometries, for example flat screens as well as more complex coil 

45 geometries, some of which are used to produce linear gradients in NMR imaging. 

On the basis of the analytical expression for the component of the magnetic field B, which is set out in 
equation (52), equations (59) and (60) become, for k- space, 
all, (ka) = -bl, (kb)j] (k) • cl> (kc)jj (k), for r>b 61 ) 
O = bK,(kb)ji (k) + cKi(kc)]J (k), for r<a (62) 

50 which yield the current densities 
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- la 1 1 (ka) (V 1 1 Oce ) K t (kb H ~ l 



(63) 



ana 

10 

jf(W - -jj(k)b Ki.(kb) 

* * ciqr^7 (64) 

For k-0, Ki{kb)/Ki(kc) = c/b and h(ka)/li(kb) = a/b from which it is deduced that the total currents It, l 2 
flowing in Si and S 2 are equaJ and opposite. From equations (63) and (64) there is obtained 

20 Mr f vZ-,-1 



F il"!, [cl I 



I, - -M-i iv-ir: i (65) 



25 = -l 2 (66) 

If a = c thereby making the primary coil and S 2 coincident, the condition b = a>/2 makes h = -I = =fe 
and means that the coil and screens may be placed in series. For discrete wire arrays chosen to 
approximate the required continuous current distributions, equations (65) and (66) become 

o o l L 1° 

(67) 

as 

and 

N,l, = -N 2 I 2 (68) 

where N 0 I 0 , Nih and N 2 I 2 are the ampere-turns for the primary coil and screens Si and S 2 respectively. A 
more general series coil arrangement is possible by varying both the turns and screen radii in equations 
4Q (67) and (68). 

A Maxwell pair of two hoops with opposite currents may be used to generate a z-gradient field. Since 
double screening produces the free space field for the hoops in the region r < a. the usual coil spacing aV3 
obtains for the most linear gradient along the z-axls. 

Figure 25 shows a simple saddle coil of radius a used to produce a transverse gradient Let this be 
45 shielded by two active cylindrical screens Si and S 2 with radii b and c respectively where a £ c < b. For a 
standard saddle geometry with 120° arcs, the primary current is 

J*(*.z) = l(5(z-dt) + «(z + di) - «(z-d 2 ) - 5(z + d 2 )} * {H(^ + »/3)[1-H(^/3)hH(^/3)[1-H(^ + 2tr/3)]} 
(69) 

where H(o) is the Heaviside function. The Fourier transform of equation (69) is 

50 

j= (k) . 23in(rrr/l) lfcgs(kd2) _ coa(kdl ,j 

0 nflT ** 

(70) 

55 

This is zero for m even or an integer of 3. Since there are now many values of m, equations (61) and (62) 
are generalised to give for two shields JJ (k)al m (ka) = -bl m (kb)j * • m (k) • cl m (kc)j J • m (k) (71 ) 
O-bK^kb)] J.™ (k) + clC (kcjj 2,m (k) (72) * • 
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From equations (71 ) and (72) there is obtained the results that 




t 



5 



co- 




rn) 



TO 



and 




00 



00 



:<l(kt) 



(74) 
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These results evaluated at k = 0 give for each separate arc pair at'tdi and td 2 the total azimuthal screen 
currents It and The dominant components of these currents arise from the m = 1 terms and may be 
25 simplified by noting that at k = 0. h'(ka)/l/(kb) = 1 and K, \kb)/Ki \kb)/Ki '(kc) = c^/b 2 . The z-components 
of current flowing in the screens may be calculated by noting that 



which expanded gives 

ajito s - i djy^0 (76) 

Fourier transforming equation (76) gives 

iTW = -frjJM (77) 

Equation (77) for R = b or c and equations (73) and (74) give on transforming to z-space the actual screen 
surface densities. 

The above results have shown that by introducing a second active screen the spatial response within a 
40 primary coil can be made to be independent of the surrounding screens. The inner screen may be 
positioned to be coincident with the primary coil and still remain the above advantages. While in the above 
description two active magnetic screens have been employed it is possible to extend the principles of 
active magnetic screening of coil structures to multiple screens. Use of two or more than two screens has 
the advantage that the screens can be designed so as not to vitiate or change the character of the magnetic 
^ field spatial response from the priamry coil structure being screened. This is true even when the inner 
screen of a two or multiple screen structure is coincident with the primary coil structure. 

The calculations and analytical expressions presented refer to continuous current density distribitions in 
the screens. Practical active screens require discrete wire arrays which simulate the continuous current 
density distributions. Discrete screens also allow exploitation of the selective transmissive and reflective 
sq properties of the active magnetic shields. 

The object of a magnetic screen is to provide a spatial current distribution which mimics that which is 
theoretically induced in a real and/or fictitious continuous conducting sleeve around a coil structure in which 
the coil itself is producing a time-dependent magnetic field gradient. Equation 24 describes an arrangement 
in which wires are equally spaced by with currents chosen again to mimic the induced surface current 
ss distribution in a continuous metal screen. 

Several methods of varying the current in these conductors are:- 

1 . To include in each conductor a small resistor chosen to give required current 

2. To change the diameter or the shape of the wire so as to affect its resistance in the right manner. 

3. Change the composition of the conductor to affect its resistance. These situations are covered in 
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Figures 30 and 31. 

It will also be clear from the above discussion that when the conductor size is changed as in Figure 31 
and in the limit of uninsulated touching wires, we have the situation shown in Figure 27. An alternative 
arrangement to this is a profiled cross-sectional band of material as shown in Figure 26. Alternative ways of 

5 producing this band or its effect are shown in Figures 28 and 29. The thickness t of the band must be 
chosen such that the electro magnetic penetration depth 5 is less than t for the highest frequency present in 
the current switching waveform. 

It will be seen that the above conductor arrangements 26 to 31 could be parallel arrangements or bands 
as in Figure 34 fed with appropriate total current along the edges indicated. However, it is conceivable that 

10 one could manufacture multi turn band structures producing a series arrangement as in Figure 35. 

An alternative approach using a constant standard wire section employs the wire arrangements in 
Figures 32 and 33. Here the wires are stacked in such a manner as to produce the desired current 
distribution. If the wires are uninsulated and touching this is an alternative method of producing the 
equivalent arrangement of Figure 26. However, if the wires are insulated and touching it will readily be seen 

is that all the turns may be in series. A series arrangement will have a much higher inductance but require 
only a small common current through each turn. The parallel arrangements discussed earlier require the 
driver circuitry to provide the total field screen current for a one turn arrangement as in Figure 34. 

In order to obtain images by magnetic resonance it is necessary to switch gradients rapidly. For typical 
magnetic gradient strengths the current I required is commonly around 150 Amps. These currents, when 

20 flowing in the static field B produce a force per unit length of F = I x B on the wires carrying this current 
For parts of the gradient coil the field B and the current I may "be perpendicular thus maximising this force. 
The resultant motion of the wires causes acoustic noise'which can be very loud if strongly coupled to the 
coil former. This problem is growing in severity with the use of higher static fields. The situation is further 
exacerbated by used of larger gradient coils where the wire length is greater over which the force can act 

25 Rapid imaging strategies can also create more noise. 

Solutions to the noise problem such as embedding the wires in rigid materials like concrete do help by 
lowering the natural resonances of the coil former and by absorbing some acoustic energy. Surrounding the 
wires with some soft acoustic absorbing material such as cotton wall can also reduce the noise. But these 
. approaches treat the symptoms rather than the cause. 

30 In a further example of the use of the present invention the problem is. solved at source by reducing the 
B field to zero in the vicinity of the wires. This eliminates the force on and hence motion of the wires. As 
well as solving the noise problem, the lack of motion of the wires removes the possibility of progressive 
fatiguing and fracture of the conductors. 

In order to do this the wires on the gradient coil are locally screened magnetically from B. This is 

35 achieved by using the principles of active screening. A single screen arrangement is created by arranging 
the screening wires in series. Different arrangements are shown in Figure 36. With these arrangements the 
far field, i.e. coordinates x. z for a point P a distance R » a where a is the half separation of the screen pair 
(Figure 36). is effectively unperturbed. For an infinite straight wire screen the interior and exterior screening 
fields B c and B p (R *> a) are respectively 



3 -a U - 

"rVa ( 78a) 



* — — 3 " > (78b) 
rrr *vr 



where R 2 = x 2 + z 2 (78c) 
*> For small separation B c can be very large and B p small. This represents a net fall off rate which goes as 
l/R 2 , For finite wires the fall off is basically dipolar i.e. goes as l/R 3 . 

For parallel infinite sheets the fields are 
B c = u^J (79a) 
Bp = 0 (79b) 

55 where J is the current density per unit length. For the arrangements sketched in Figures 36a to 36c 
respectively the screening efficacy will lie between the cases covered in Equations. (78) and (79). Because 
the screen produces a static field, perturbation of the main field can be eliminated with a shim set The 
screen or counter field generator need be active only during the experimental period. 
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To further reduce the extraneous static fields produced by the wire screening arrangement, a second 
active screen may be used in addition to the first screen. 

Such a double active screen arrangement is sketched in Figure 36. The wire W is roughly screened by 
a parallel plate arrangement Si Residual far fields . arising from non-cancellation outside the plates is 
5 annulled by a second active screen S2 comprising a set of conductors distributed appropriately. Figure 37 
shows two local wire screening strategies for a circular hoop (37a) and a saddle gradient (37b). 

Using the wife screening arrangement described, the benefits of rapid switching of large gradients may 
be obtained within a high static magnetic field without acoustic noise. 

Screening could also be applied to a HaJI probe to increase its sensitivity to small field variation by 
to removing the main central field. 

In addition to the application of active magnetic screening to electromagnet and gradient coil design 
there are a number of possible applications in rf coil design. For example a fully screened series wound rf 
coil placed coaxial ly inside an unscreened coil would have zero mutual inductance. However an NMR 
sample placed inside the screen coil would sense Melds generated by both coils. The coils though coaxial 
15 therefore behave electrically as though they were orthogonal. This can have advantages in multi-nuclear 
irradiation and detection. 

A schematic example of one such possible configuration is shown in Figure 38 where an inner coil C1 
is wound inside an outer coil C2. The inner coil C1 can be screened by providing a suitably wound 
screening coil CS between the two coils C1 and C2. The screen coil CS is as shown connected in series 
20 with coil C1 to pass the same current It is generally longer than the coils C1 and C2 to provide efficient 
screening and is wound such that the current in coil CS opposes that in coil C1. The exact positioning of 
the wires in coil CS is determined by using the above approach. 



25 Claims 

1. A gradient coil set for an MR apparatus comprising a plurality of coils, each of said coils adapted to 
provide a respective surface current distribution, the magnetic field resulting from said surface current 
distributions exhibiting a predetermined area inside said coil set and a substantially zero value outside said 

30 coil set 

2. The coil set of claim 1 wherein each of said coils is cylindrically-shaped. 

3. The coil set of claim 1 wherein said coils are connected in series, said coil set being adapted to 
receive its operating current from a single power supply. 

4. A coil system for producing a magnetic field comprising: 

35 a first coil adapted to be coupled to a source of current for providing a first surface current on the surface of 
a first cyliner having a radius n and an axial length Zi ; and 

a second coil adapted to be coupled to a source of current for providing a second surface current on the 
surface of a second cylinder substantially coaxial with said First cylinder and having a radius r 0 greater than 
n and an axial length zo; 

40 said first and second surface currents mutually providing a magnetic field having a predetermined gradient 
in a predetermined dimension and a substantially constant value across the other two dimensions within the 
volume enclosed by said first coil and having a substantially zero value in the volume outside of said 
second coil. 

5. The coil system recited in claim 5 wherein said axial length zo is greater than said axial length zi . 
45 6. The coil system of claim 4 wherein said coils are connected in series. 

7. A magnetic resonance system comprising: 
a main magnet for establishing a uniform magnetic field in an imaging volume in the interior of said main 
magnet 

radio-frequency means for pulsing said imaging volume with electromagnetic energy for stimulating nuclear 

so magnetic resonance of nuclei within said imaging volume; 

detecting means for sensing the nuclear magnetic resonance of nuclei within said imaging volume; and 
a first gradient coil set located within said uniform magnetic field of said main magnet, said first coil set 
being comprised of an inner and an outer coil, said inner coil adapted to be coupled to a source of current 
for providing a first surface current distribution on the surface of a first cylinder having a radius n an an 

55 axial length zi, said outer coil adapted to be coupled to a source of current for providing a second surface 
current distribution on the surface of a second cylinder substantially coaxial with said cylinder and having a 
radius r 0 greater than n and an axial length zo, said first and second surface currents being adapted to 
produce a magnetic field having a gradient in a first dimension and a constant value across the other two 
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dimensions within said first coil set and having a substantially zero value outside of said first coil set. 

8. The magnetic resonance system of claim 1 1 further comprising second and third gradient coil sets 
for establishing magnetic fields having gradients in each of the two remaining dimensions, respectively, 
within the respective cot! sets and having substantially zero values outside the respective coil sets. 
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